REPORT  82-39 

ADA  124516 


US  Army  Corps 
of  Engineers 

Cold  Regions  Research  & 
Engineering  Laboratory 


The  brine  zone  in  the  McMurdo  Ice  Shelf, 
Antarctica 


For  conversion  of  SI  metric  units  to  U.S.I 
British  customary  units  of  measurement 
consult  ASTM  Standard  E380,  Metric  Prac¬ 
tice  Guide,  published  by  the  American  Socie¬ 
ty  for  Testing  and  Materials,  1916  Race  St., 
Philadelphia,  Pa.  19103. 


Cover:  Coring  at  site  of  brine  step  on  the 
McMurdo  Ice  Shelf  in  January  1977. 
Mt.  Erebus  is  seen  in  the  background. 


CRREL  Report  82-39 

December  1982 


The  brine  zone  in  the  McMurdo  Ice  Shelf, 
Antarctica 

A.  Kovacs,  A.J.  Gow,  J.H.  Cragin  and  R.M.  Morey 


Prepared  for 

NATIONAL  SCIENCE  FOUNDATION 
Approved  for  public  release;  distribution  unlimited. 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whan  Data  Entarad) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.  REPORT  NUMBER  2.  GOVT  ACCESSION  NO. 

CRREL  Report  82-39  H i/H  5/L 

3.  RECIPIENT'S  CATALOG  NUMBER 

9 

4.  TITLE  (and  Subtitle) 

THE  BRINE  ZONE  IN  THE  McMURDO 

ICE  SHELF,  ANTARCTICA 

5.  TYPE  OF  REPORT  A  PERIOD  COVERED 

6.  PERFORMING  ORG.  REPORT  NUMBER 

7.  author^*; 

A.  Kovacs,  A.J.  Gow,  J.H.  Cragjn  and  R.M.  Morey 

8.  CONTRACT  OR  GRANT  NUMBERfcj 

NSF-DPP  74-23654 

NSF-DPP  8004221 

*-  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

U.S.  Army  Cold  Regions  Research  and  Engineering  Laboratory 

Hanover,  New  Hampshire  03755 

10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  ft  WORK  UNIT  NUMBERS 

•  1.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Division  of  Polar  Programs 

National  Science  Foundation 

Washington,  D.C.  20314 

12.  REPORT  DATE 

December  1982 

13.  NUMBER  OF  PAGES 

35 

14.  MONITORING  AGENCY  NAME  4  AODRESSffl  dll  (a  rant  tram  Controlling  Olllea) 

IS.  SECURITY  CLASS,  (ol  thla  roport) 

Unclassified 

15a.  declassification/downgrading 
SCHEDULE 

IS.  DISTRIBUTION  STATEMENT  (at  thla  Report) 

Approved  for  public  release;  distribution  unlimited. 

17.  DISTRIBUTION  STATEMENT  (ol  lha  abatract  an  farad  In  Block  30,  II  dlllaranl  from  Roport) 

IS.  SUPPLEMENTARY  NOTES 

IS.  KEY  WORDS  (Contlnua  on  tararaa  alda  II  nocoooary  and  Identity  by  block  number) 

Antarctica  Ice  shelves 

Brines  Radio  echo  sounding 

Firn 

Glaciers 

2G  ABSTRACT  (Centime  am  rararaa  at*  It  mnnaqr  and  Identity  by  block  mmtbor) 

N  Observations  of  a  4.4-m-high  brine  step  in  the  McMurdo  Ice  Shelf,  Antarctica,  show  that  it  has  migrated  about  1 .2 
km  in  4  years.  The  present  brine  wave  is  overriding  an  older  brine-soaked  layer.  This  migration  is  proof  of  the  dynamic 
nature  of  the  step,  which  is  the  leading  edge  of  a  brine  wave  that  originated  at  the  shelf  edge  after  a  major  break-out  of 
the  McMurdo  Ice  Shelf.  The  inland  boundary  of  brine  penetration  is  characterized  by  a  series  of  descending  steps  that 
are  believed  to  represent  terminal  positions  of  separate  intrusions  of  brine  of  similar  origin.  The  inland  boundary  of 
brine  percolation  is  probably  controlled  largely  by  the  depth  at  which  brine  encounters  the  fim/ice  transition  (43  m). 
However,  this  boundary  is  not  fixed  by  permeability  considerations  alone,  since  measurable  movement  of  brine  is  still 
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THE  BRINE  ZONE  IN  THE 
MCMURDO  ICE  SHELF,  ANTARCTICA 


A.  Kovacs,  A.J.  Gow,  J.H.  Cragin  and  R.M.  Morey 


INTRODUCTION 

The  occurrence  of  brine-soaked  firn  in  Antarc¬ 
tic  ice  shelves  has  been  reported  by  Dubrovin 
(1962)  for  the  Lazarev  Ice  Shelf,  Stuart  and  Bull 
(1963)  for  a  small  ice  shelf  near  Wilkes  Station, 
Smith  and  Evans  (1972)  for  the  Wordie,  Brunt  and 
Larsen  Ice  Shelves,  and  Thomas  (1973)  for  the 
Brunt  Ice  Shelf.  This  phenomenon,  however,  is 
best  known  for  the  McMurdo  Ice  Shelf  as  report¬ 
ed,  for  example,  by  Stuart  and  Bull  (1963),  Risk 
and  Hochstein  (1967),  Heine  (1968),  Clough 
(1973),  and  Kovacs  and  Gow  (1975,  1977b).  Sev¬ 
eral  mechanisms  have  been  suggested  for  the  oc¬ 
currence  of  brine  soaking,  including  vertical  per¬ 
colation  of  seawater  from  the  bottom,  lateral  infil¬ 
tration  from  the  shelf  edge,  and  upward  diffusion 
along  ice  crystal  boundaries.  The  most  recent 
work  by  Kovacs  and  Gow  (1975,  1977b)  tends  to 
favor  a  process  of  lateral  infiltration  into  dry  firn 
by  seawater.  Essential  to  this  concept  of  infiltra¬ 
tion  of  ice  shelves  by  seawater  is  the  need  for  sub¬ 
stantially  rapid  elimination  of  excess  water  by 
freezing  within  pore  spaces  in  the  permeable  firn. 
This  process  leads  to  the  formation  of  concen¬ 
trated  brine  capable  of  remaining  liquid  at  the 
temperatures  encountered  in  the  infiltrated  firn. 
Such  a  process  is  also  accompanied  by  changes  in 
brine  chemistry  (Wilson  and  Heine  1964,  Heine 
1968,  Cragin  et  al.  1982)  and  very  substantial  iso¬ 
topic  exchange  between  the  brine  and  the  firn 
(Stewart  1975). 


OBJECTIVES 

The  major  objectives  of  this  study  were  to  delin¬ 
eate  accurately  the  terminal  (inland)  boundary  of 
seawater  infiltration  in  the  McMurdo  Ice  Shelf 
and  to  examine  the  physical  structure,  chemistry, 
depth  characteristics  and  lateral  continuity  of  the 
brine  layer.  These  investigations  combined  radio 
echo  profiling  with  core  drilling  at  sites  selected  on 
the  basis  of  the  radio  profiling  results.  Coring  in 
this  manner  permitted  direct  examination  of  the 
brine  infiltrated  and  noninfiltrated  ice  shelf  struc¬ 
ture.  Studies  conducted  on  the  cores  included 
chemical  analysis  of  brine-impregnated  samples  as 
a  means  of  monitoring  the  freeze-concentration 
of  the  seawater  as  it  permeated  the  ice  shelf  firn, 
density  profile  measurements,  and  petrographic 
examinations  of  the  structure  of  the  brine-infil¬ 
trated  firn.  The  ultimate  objective  was  to  resolve 
four  questions:  1)  What  are  the  principal  mechan¬ 
isms  by  which  seawater  infiltrates  the  McMurdo 
Ice  Shelf?  2)  What  are  the  structural  relationships 
of  the  brine  zone  to  those  parts  of  the  ice  shelf  that 
remain  unaffected  by  the  infiltration  of  seawater? 
3)  Why  does  the  infiltration  of  brine  terminate 
where  it  does?  4)Is  this  termination  of  brine  mi¬ 
gration  mainly  controlled  by  permeability  consid¬ 
erations  or  are  other  factors  involved? 

The  study  was  confined  to  the  portion  of  the 
Ross  Ice  Shelf  generally  designated  the  McMurdo 
Ice  Shelf  (see  Fig.  I).  The  part  affected  by  brine 
infiltration  is  included  within  the  area  of  net  sur- 
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Figure  I.  Sketch  map  of  McMurdo  Sound  and  environs 
showing  the  portion  of  the  Ross  Ice  Shelf  designated  the 
McMurdo  Ice  Shelf  (dotted  area). 


face  snow  accumulation  and  shelf  bottom  ablation 
which  occurs  mainly  in  the  region  directly  east  and 
south  of  Hut  Point  Peninsula.  Further  west  the  ice 
shelf  regime  gives  way  to  net  surface  ablation  and 
bottom  freezing. 


ANALYTICAL  TECHNIQUES 
Radio  echo  profiling 

The  depth  characteristics,  lateral  continuity  and 
inland  boundary  of  seawater  infiltration  into  the 


McMurdo  Ice  Shelf  were  determined  in  January 
1977  with  a  dual-antenna  radio  echo  profiling  sys¬ 
tem.  The  signal  information,  i.e.  the  two-way 
travel  time  of  the  transmitted  wavelet  to  and  from 
various  reflecting  interfaces,  was  displayed  on  a 
graphic  recorder  and  stored  on  a  magnetic  tape 
recorder  for  later  playback  and  analysis.  The  dual¬ 
antenna  arrangement  is  depicted  in  Figure  2  and 
an  example  of  a  received  waveform  is  shown  in 
Figure  3.  With  the  dual-antenna  arrangement,  it 
was  possible  to  determine  the  velocity  of  propaga¬ 
tion  of  the  electromagnetic  (EM)  wave  in  the  firn 
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Figure  2.  Schematic  of  dual-antenna  configura¬ 
tion  and  transmitted  wa-  '?/  paths. 
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Figure  3.  Example  of  received  radio  echo  wavelet.  In  the 
graphic  record  many  horizons  or  internal  layers  were  observed  in  the 
ice  shelf  firn.  These  layers  could  be  tracked  for  many  kilometers. 
Their  voltage  amplitude  signature  is  shown  above. 
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Figure  4.  Dual  antenna  radio  echo  profiling  arrangement  (a)  and  snow  pit 
studies  at  site  of  4.4-m-high  brine  step  at  station  10  (b). 


and  the  firn’s  effective  bulk  dielectric  constant, 
irrespective  of  the  firn’s  changing  properties  due 
to  density. 

One  antenna  operated  in  a  transmit-receive 
mode  and  the  second  in  a  receive-only  mode.  The 
antennas  were  placed  a  fixed  distance  apart  on  the 
snow  surface  and  moved  as  a  unit  (see  Fig.  4).  The 
effective  velocity  ( Ke)  of  the  EM  wave  was  then  de¬ 
termined  by: 


K  = 


0) 


where  x  = 
'd  = 


t 


a 


distance  between  centers  of  the  two 
antennas 

vertical  travel  time  from  transceiver 
antenna  to  and  from  the  subsurface 
interface 

travel  time  from  transceiver  antenna 
to  subsurface  interface  to  receive- 
only  antenna  (rx .)  plus  time  ra 
air-wave  travel  time  between  anten¬ 
nas  (see  Fig.  2). 


A  more  detailed  description  of  the  methodology 
associated  with  use  of  the  dual-antenna  sounding 
system  can  be  found  in  Kovacs  and  Morey  (1979). 

With  the  effective  velocity  of  propagation  thus 
determined,  it  is  then  possible  to  determine  brine 
layer  depth  D  from 


and  the  real  effective  bulk  dielectric  constant  ere  of 
the  firn  from 


'df 

2D 


(3) 


where  c  is  the  free  space  electromagnetic  velocity. 

In  addition  to  furnishing  a  continuous  trace  of 
the  top  of  the  brine  layer  the  radio  echo  profiler 
can  also  detect  such  features  as  internal  firn  layer¬ 
ing,  cracks,  and  relict  brine  horizons  (Kovacs  and 
Gow  1975,  1977a).  However,  it  should  be  empha¬ 
sized  here  that  the  major  discontinuity  observed  in 
radio  echo  profiling  records  in  the  McMurdo  Ice 
Shelf  is  due  entirely  to  a  signal  return  from  the  top 
of  the  brine  layer.  This  situation  has  been  inde¬ 
pendently  confirmed  by  drilling  and  coring  at 
seven  locations  on  the  McMurdo  Ice  Shelf.  In  all 
our  drill  holes  we  never  encountered  rock  debris  or 
any  other  englacial  material  capable  of  creating 
such  a  strong  reflecting  horizon  with  the  ice  shelf. 


Core  drilling 

Drill  holes  were  all  located  along  the  line  of  the 
main  radio  echo  profiling  traverse  of  January 
1977.  An  electromechanical  drill  (Fig.  5),  designed 
by  Rand  (1976)  and  operated  by  a  team  of  drillers 
from  the  Polar  Ice  Coring  Office  (PICO),  was 
used  for  deep  coring  at  four  sites.  Three  of  these 
holes  penetrated  the  brine  layer  for  distances  of  up 
to  2.5  m,  but  none  were  drilled  deep  enough  to 
penetrate  the  entire  brine  zone.  The  fourth  and 
deepest  hole  was  purposely  drilled  a  few  meters 
beyond  the  brine  terminus,  as  determined  by  radio 
echo  profiling,  so  as  not  to  intercept  the  brine 
zone.  Cores  from  all  four  holes  were  logged  and 
measured  on-site  for  temperature.  Those  cores  in¬ 
filtrated  with  brine  were  cut  into  =  10-cm  lengths 
and  placed  in  sealed  containers.  Cores  and  sam¬ 
ples  were  then  transported  back  to  refrigerated 
facilities  at  McMurdo  Station,  where  the  brine- 
soaked  samples  were  first  measured  and  weighed 
for  density  determination  and  then  melted.  The 
salinity  of  the  meltwater  from  each  sample  was 
then  measured  and  specimens  of  the  liquid  bottled 


Figure  5.  Electromechanical  drill  hoist  assembly 
used  for  coring  into  the  ice  shelf. 
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RESULTS  AND  DISCUSSION 


and  carefully  sealed  for  subsequent  chemical  anal¬ 
ysis  at  CRREL.  Thin  sections  of  selected  samples 
were  also  prepared  and  examined  cursorily  for 
crystal/air  bubble/brine  inclusion  relationships. 
Unfortunately  all  these  sections,  together  with  the 
remaining  brina-soaked  samples  set  aside  for  addi¬ 
tional  structural  studies,  were  lost  when  the 
McMurdo  Station  coldroom  in  which  they  were 
stored  was  accidentally  defrosted.  Cores  from 
several  shallow  holes  located  near  the  front  of  the 
McMurdo  Ice  Shelf  were  also  investigated  in  con¬ 
junction  with  cores  from  the  deeper  holes.  Unused 
cores  were  returned  to  CRREL  for  archiving  and 
additional  analysis. 


Brine  infiltration  survey 
During  January  1977  radar  profiling  of  the  top 
of  the  brine  layer  was  conducted  along  several 
separate  traverses,  including  the  main  profile  line, 
which  was  flagged  and  resurveyed  in  November 
1978  and  January  1981.  The  main  profile  line  ex¬ 
tended  from  the  ice  shelf  edge  to  the  inland 
boundary  of  brine  infiltration  and  for  some 
distance  beyond.  In  addition,  the  inland  boundary 
of  brine  infiltration  was  profiled  in  detail  between 
Hut  Point  Peninsula  and  White  Island  (Fig.  6). 
Stations  along  the  boundary  were  fixed  by  tri- 


Figure  6.  Area  encompassed  by  radio  echo  profiling  surveys.  The  route  of  the  major 
profiting  traverse  and  the  boundary  marking  the  inland  limit  of  brine  infiltration  are  indi¬ 
cated,  together  with  pertinent  data  on  ice  shelf  movement  and  surface  features.  The  boun¬ 
dary  marking  the  westward  limit  of  brine  infiltration  (the  equilibrium  line)  is  also  indi¬ 
cated.  The  shelf  to  the  west  of  this  boundary  is  ablating  glacial  ice.  The  area  inside  the  first 
dashed  line  immediately  west  of  White  Island  is  composed  of  snow  ice  up  to  15  m  thick. 
Snow  depth  is  less  than  2.5  nt. 


Table  1. 
ice-shelf 

Station 

no. 

1977  McMurdo  Ice  Shelf  station  data  along  traverse  from  ice  front  to  beyond 
movement  marker  307. 

Effective 

Keel  EM  wave 

Distance  Freeboard  Brine  Brine  depth  Real  effective  velocity 

inland  F  depth  elevation  K  dielectric  constant 

On)  On)  On)  f-tn)  (in)  <rp  (m/ns) 

i 

40 

3.9 

4.7 

0.8 

1.81** 

0.223 

2 

305 

4.8 

— 

— 

15.5* 

— 

— 

3 

600 

5.0 

6.0 

1.0 

— 

1.91 

0.217 

4 

750 

5.6 

— 

— 

20.0* 

— 

— 

5 

910 

6.0 

— 

— 

19.5* 

— 

— 

6 

1.080 

6.5 

7.8 

1.3 

21.5* 

1.95 

0.225 

7 

1.640 

8.1 

10.2 

2.1 

— 

1.99 

0.212 

8 

2.160 

10.2 

12.8 

2.6 

— 

2.04 

0.210 

9 

2,620 

11.6 

15.5 

3.9 

— 

2.09 

0.208 

10 

2.840 

12.5 

16.9 

4.4 

— 

2.10 

0.207 

10b 

2.900 

12.7 

21.5 

8.8 

— 

2.18 

0.203 

II 

3,280 

14.4 

23.6 

9.2 

— 

2.18 

0.203 

12 

3,940 

15.6 

25.4 

9.8 

— 

2.22 

0.201 

13 

4,550 

16.7 

26.7 

10.0 

— 

2.25 

0.200 

14 

5,190 

17.6 

27.9 

10.3 

— 

2.29 

0.198 

15 

5,640 

18.2 

29.2 

11.0 

— 

2.32 

0.197 

16 

6,160 

19.4 

31.2 

11.8 

— 

2.33 

0.196 

17 

6,680 

20.1 

32.6 

12.5 

— 

2.34 

0.196 

18 

7,200 

20.8 

33.8 

13.0 

— 

2.35 

0.196 

19 

7,720 

21.7 

35.3 

13.6 

— 

2.36 

0.195 

20 

8,480 

23.0 

37.9 

14.9 

— 

2.43 

0.192 

21 

high 

8,800 

23.4 

40.6 

17.2 

— 

2.46 

0.191 

21 

low 

8,800 

23.4 

44.5 

21.1 

— 

2.50 

0.190 

22 

high 

9,080 

23.8 

45.5 

21.7 

— 

2.49 

0.190 

22 

low 

9,080 

23.8 

47.5 

23.7 

— 

2.52 

0.198 

•  Drill-hole  measurements. 

t  Based  on  wavelet  two-way  travel  time  from  snow  surface  to  shelf  bottom  minus  two-way  travel  time  from 
snow  surface  to  brine  elevation  at  station  24.  The  ice  below  brine  termination  elevation  was  assumed  to  have 
f  =  2.95  and  therefore  Vt  =  0.175  m/ns. 

•*  Between  surface  and  brine  layer. 

Note:  Brine  elevation  is  distance  from  sea  level  to  top  of  brine  layer.  Keel  depth  is  draft  of  ice  shelf. 


angulation  survey.  The  resulting  data  were  com¬ 
puter-analyzed  for  accurate  plotting  on  the  New 
Zealand  survey  grid  of  the  area.  The  northern 
limit  of  this  brine  layer  was  originally  delineated 
by  Clough  and  Bentley  (1967)  and  Clough  (1973), 
who  emplaced  bamboo  pole  markers  at  intervals 
along  the  boundary.  In  November  1973,  new  bam¬ 
boo  poles  were  placed  at  intervals  along  the  boun¬ 
dary  (John  Clough,  pers.  comm.).  The  18  January 
1977  resurvey  of  the  existing  markers  2S,  IS,  IN, 
2N,  3N  and  4N  (Fig.  6)  indicated  that  the  brine 
terminus  had  moved  inland  197,  129,  73,  S3,  10 
and  9  m  beyond  these  stations  respectively  during 


the  intervening  3  years,  i.e.  at  average  rates  of 
0.168  ,  0.110,  0.062,  0.047,  0.009  and  0.008 
m/day.  The  magnitude  of  brine  migration  at  the 
first  four  markers  demonstrates  that  brine  infiltra¬ 
tion  there  is  still  very  active.  However,  the  very 
small  movement  recorded  at  markers  3N  and  4N 
would  indicate  that  the  brine  there  is  now  in¬ 
filtrating  firn  with  very  limited  permeability. 

The  brine  layer  depth  and  the  effective  velocity 
of  propagation  of  the  EM  wave  and  real  effective 
dielectric  constant  values  for  the  firn,  as  deter¬ 
mined  from  the  ratio  echo  profile  data,  versus  dis¬ 
tance  along  the  main  profile  line  shown  in  Figure  6 
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are  listed  in  Table  1.  Also  listed  are  the  surface 
elevations  (freeboards)  determined  by  topographic 
survey  and  the  ice  shelf  keel  depths  determined  by 
direct  drill  hole  measurement.  These  data  allowed 
the  construction  of  a  representative  cross  section 
of  the  McMurdo  Ice  Shelf  along  our  main  profile 
line  (Fig.  7). 

Brine  layer  steps 

Among  several  important  features  shown  in 
Figure  7  are  a  large  4.4-m  step  in  the  brine  layer  lo¬ 
cated  near  station  10  and  a  series  of  descending 
steps  characterizing  the  inland  termination  of  the 
brine  layer.  These  same  features,  as  they  appear 
on  the  graphic  recorder,  are  illustrated  in  Figures  8 
and  9. 

The  position  of  the  4.4-m  step  in  the  brine  layer, 
as  first  observed  on  18  January  1977,  was  careful¬ 
ly  flagged  and  reprofiled  on  24  November  1978 
(Fig.  7)  and  again  on  20  January  1981.  Between  18 
January  1977  and  24  November  1978  this  step  was 
found  to  have  moved  inland  about  700  m,  equiva¬ 
lent  to  an  average  rate  of  advance  of  about  1 .04 
m/day.  Between  24  November  1978  and  20  Janu¬ 
ary  1981  the  brine  step  had  migrated  an  additional 
530  m  inland  at  an  average  rate  of  0.67  m/day. 
These  measurements  of  the  actual  migration  of  the 
brine  step  demonstrated  for  the  first  time  the  dy¬ 
namic  wave-like  nature  of  seawater  infiltration 
within  the  McMurdo  Ice  Shelf.  This  particular 
brine  step  almost  certainly  represents  the  leading 
edge  of  a  wave  of  seawater  that  entered  at  the  ice 
front  during  a  recent  major  break-out  of  ice  in 
McMurdo  Sound.  In  addition,  it  would  appear 
that  over  the  4-year  duration  of  the  measure¬ 
ments,  brine  infiltration  has  slowed  appreciably, 
most  probably  in  response  to  a  decrease  in  firn 


permeability  as  the  brine  moved  into  deeper  and 
denser  parts  of  the  ice  shelf. 

Results  of  drilling  at  stations  10  and  10b  on 
either  side  of  the  brine  step  clearly  show  that  the 
4.4-m-high  brine  wave  is  currently  riding  on  top  of 
an  earlier  brine-soaked  layer.  This  layer  also  ter¬ 
minates  in  a  step,  some  distance  inland,  and  it,  to¬ 
gether  with  several  other  steps  located  near  the 
present  boundary  of  brine  infiltration,  are  also  be¬ 
lieved  to  represent  earlier  episodes  of  wave-like  in¬ 
trusions  of  seawater  originating  at  the  ice  front. 
All  such  intrusions  of  seawater  are  believed  to  be 
triggered  by  periodic  break-outs  of  the  ice  shelf. 
Through  this  process  firn  at  the  shelf  edge  in  a 
zone  between  sea  level  and  the  top  of  the  old  brine 
horizon,  now  located  below  sea  level,  becomes  ex¬ 
posed  to  the  sea  and  a  new  wave  of  seawater  be¬ 
gins  to  percolate  into  the  ice  shelf.  This  sequence 
of  events  is  illustrated  in  Figure  10.  From  the 
measured  slope  of  the  brine-soaked  firn  layer,  it 
was  estimated  that  approximately  3  km  of  the  ice 
shelf  had  to  have  broken  off  to  allow  for  the  for¬ 
mation  of  the  4.4-m-high  brine  step.  Paige  (1971) 
reports  that  “the  most  extensive  observed  periods 
of  ice  shelf  breakout  occurred  in  February  of  1964 
and  1965.”  Approximately  1500  m  of  shelf  calved 
during  these  events.  Further  calving  occurred  in 
1966,  1967,  1968,  and  apparently  in  1970  but  these 
were  less  extensive.  Nevertheless,  the  sum  total  of 
the  1964  through  1970  break-outs  must  have  been 
about  3  km.  Since  no  brine  step  record  of  these 
events  was  observed  in  the  radio  echo  profiles  we 
can  only  assume  that  these  steps  were  removed  by 
subsequent  break-outs.  If  ice  shelf  break-out  oc¬ 
curred  in  February  1970  and  the  resultant  brine 
wave  coincided  with  the  4.4-m  step  located  2850  m 
in  from  the  shelf  edge  as  of  January  1977,  then  the 


Figure  8.  Graphic  record  of  1977  radio  echo  data  showing  the 
4.4-m-high  brine  step  located  at  station  10. 
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a.  Time  between  horizontal  scan  lines  is  about  325  ns  and  the  distance  between  station  E  and  F  is 
=  200  m. 
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b.  12.0  km  of  the  ice  shelf  bottom  is  shown,  which  reveals  a  relatively  smoothly  sloping  surface. 


Figure  9.  Graphic  record  of  1977  radio  echo  data  showing  the  stepped  nature  of  several  brine  layers 
located  near  the  inland  boundary  of  brine  penetration.  The  fact  that  the  shelf  bottom  could  be  observed  under 
the  lowest  brine  layer  indicates  that  the  layer  is  "thin”  ( 2-3  m  or  less)  and/or  low  loss  with  respect  to  electromag¬ 
netic  energy  transmission;  otherwise  the  bottom  radio  echo  would  not  have  been  observed.  . 


Before  Break-out  After  Breok-out 


Figure  10.  Schematic  of  ice  shelf  break-out  event  leading  to  formation  of  new  zone  of 
brine-infiltrated  firn. 


average  infiltration  velocity  of  the  brine  wave 
would  have  been  about  1.2  m/day.  This  rate  of  in¬ 
filtration  is  compatible  with  the  value  of  1.04 
m/day  obtained  by  direct  measurement  of  the 
4.4-m  step  between  18  January  1977  and  24 
November  1978. 

It  was  learned  from  a  discussion  with  Mr.  Cass 
Roper,  who  was  at  the  New  Zealand  Scott  Base  at 
the  time,  that  between  6  February  and  22  March 
1980,  large  sections  of  the  McMurdo  Ice  Shelf 
broke  off.  During  one  period,  the  shelf  was  ob¬ 
served  moving  vertically  3  to  5  cm  at  an  inland 
crack.  A  resurvey  of  the  main  profile  line  in  Janu¬ 
ary  1981  revealed  that  1640  m  of  the  shelf  was  then 
missing.  Based  upon  old  records  (Heine  1963, 
Paige  1971)  this  break-out  was  one  of  the  largest 
ever  recorded.  A  new  brine  layer  horizon  was  de¬ 
tected,  situated  approximately  1.4  m  above  the  old 
horizon  at  the  shelf  edge  and  extending  inland  375 
m,  as  shown  in  Figure  7.  This  infiltration  occurred 
over  a  period  of  about  302  days,  equivalent  to  an 
average  rate  of  brine  migration  of  1.24  m/day. 

Brine  infiltration  characteristics 

Risk  and  Hochstein  (1967)  attribute  the  forma¬ 
tion  of  brine-soaked  firn  in  the  McMurdo  Ice 
Shelf  to  lateral  percolation  of  seawater  from  the 
shelf  edge.  They  also  embrace  the  idea  of  Stuart 
and  Bull  (1963),  who  suggest  that  vertical  brine  in¬ 
filtration  occurs  through  low  density  bottom  firn. 
However,  Risk  and  Hochstein  limited  the  area 
where  the  latter  mechanism  occurs  to  a  zone  ex¬ 
tending  about  0.5  km  in  from  the  shelf  edge  where 
permeable  firn  was  considered  to  be  exposed  due 
to  bottom  melting.  If  vertical  infiltration  is  occur¬ 
ring  in  a  zone  near  the  shelf  edge,  one  would  ex¬ 
pect  the  top  of  the  brine-soaked  firn  in  this  area  to 
be  in  near-hydrostatic  equilibrium  with  the  sea¬ 
water.  From  our  radio  echo  profile  surveys,  it  was 
observed  that  the  top  of  the  brine  layer,  beginning 
at  the  very  edge  of  the  ice  shelf,  slopes  down  below 
sea  level,  that  is,  has  a  negative  gradient.  It  there¬ 
fore  appears  that  permeable  firn  does  not  now  ex¬ 
ist  on  the  bottom  of  the  ice  shelf.  In-situ  freezing 
of  water  derived  from  the  brine  as  it  laterally  in¬ 
filtrates  the  firn  has  effectively  blocked  off  inter¬ 
connected  pores,  preventing  any  upward  migra¬ 
tion  of  seawater. 

It  is  interesting  to  note  that  Thomas  (1975)  esti¬ 
mated  a  brine  infiltration  rate  of  1.10  m/day  for 
the  Brunt  Ice  Shelf,  where  the  firn  density  (0.570 
Mg/mJ),  temperature  (-  10°C)  and  pressure  gra¬ 
dient  (1  in  525)  were  more  favorable  to  seawater 
infiltration  than  for  the  McMurdo  Ice  Shelf.  We 
can  calculate  a  brine  flow  velocity  u  at  station  10 


following  the  procedure  of  Thomas,  where 


and  B0  is  the  specific  permeability  of  the  firn,  ij  is 
the  dynamic  viscosity  of  the  brine  fluid,  and 
dp/dx  is  the  pressure  gradient  of  the  brine  layer. 
Thomas  modified  Shimizu’s  (1970)  formula  B0  = 
0.077  xd1  exp(  -  7.4 gs),  where  d  is  the  grain  diam¬ 
eter  and  q5  is  the  density  of  snow,  and  has  sug¬ 
gested  that  B0  =  9xl0'4  d2.  We  estimated  the 
grain  diameter  to  be  2  mm  for  0.7-Mg/m’  density 
firn  at  station  10.  The  firn  temperature  was  * 
-  14  °C  and  tj  for  the  brine  at  this  temperature  is 
*  4.5  x  10"’  Pa/s.  The  hydraulic  head,  over  a  dis¬ 
tance  of  2580  m  from  station  1  to  station  9,  is  3.1 
m.  From  eq  4  we  calculate  a  flow  velocity  of 
9.68x10“*  m/s  or  *  0.84  m/day.  This  value  fits 
nicely  between  the  18  January  1977-24  November 
1978  average  brine  velocity  of  1.04  m/day  which 
existed  before  the  brine  wave  reached  station  10, 
and  the  24  November  1978-20  January  1981  aver¬ 
age  velocity  of  0.67  m/day  measured  after  the 
wave  passed  station  10. 

Core  data  obtained  in  November  1978  from 
holes  drilled  at  stations  E  and  F,  on  either  side  of 
the  brine  terminus,  revealed  that  “impermeable” 
firn  of  density  0.82-0.83  Mg/m’  (ice  by  definition) 
is  encountered  at  a  depth  of  about  43  m,  or  6  to  7 
m  above  the  brine  layer  at  the  station  E  hole, 
which  was  located  about  120  m  inside  the  brine 
terminus  as  measured  in  January  1977.  The  hole  at 
station  F  was  located  about  77  m  beyond  the  brine 
terminus,  as  it  existed  in  January  1977,  and  no 
brine  was  encountered  in  this  59-m-deep  hole 
when  it  was  drilled  in  November  1978.  From  these 
observations,  it  might  appear  that  the  inland  limit 
of  brine  infiltration  in  the  McMurdo  Ice  Shelf  is 
controlled  mainly  by  the  depth  at  which  brine  en¬ 
counters  the  firn/ice  transition.  At  this  boundary 
ice  containing  liquid  brine  inclusions  is  presum¬ 
ably  carried  downward  into  the  ice  shelf  and  fur¬ 
ther  densified  as  snow  continues  to  accumulate  on 
the  ice  shelf  surface. 

Some  estimate  of  the  time  that  has  elapsed  since 
the  brine  became  “locked  in”  can  be  obtained 
from  accumulation  data  reported  by  Heine  (1967), 
which  show  that  snow  accumulation  in  the  vicinity 
of  the  brine  terminus  should  average  around  25 
cm  water  annually.  At  station  E,  where  the  top  of 
the  brine  layer  is  now  located  6  to  7  m  below  the 
firn/ice  transition,  we  calculate,  on  the  basis  of 
our  core  density  measurements,  that  the  infiltra¬ 
tion  of  brine  should  have  effectively  ceased  20  to 


25  years  ago.  However,  more  recent,  very  careful 
profiling  of  this  boundary  has  revealed  that  a 
small  but  detectable  migration  of  brine  is  still  oc¬ 
curring,  apparently  by  some  mechanisms  not  di¬ 
rectly  related  to  firn  permeability.  Once  “locked 
in,”  the  brine  must,  with  time,  descend  into 
deeper  and  warmer  ice.  In  doing  so,  the  concen¬ 
trated  brine  would  need  to  dissolve  ice  to  maintain 
a  brine  solution  concentration  in  equilibrium  with 
the  in-situ  temperature,  and  this  may  have  con¬ 
tributed  to  the  continued  slow  inland  percolation 
of  the  brine. 

The  slope  of  the  brine  layer  measured  1  in  832 
from  the  shelf  edge  to  station  9.  From  station  1 1 
to  station  20  the  brine  layer  slope  was  1  in  910,  and 
from  stations  22  to  24  it  had  decreased  to  approxi¬ 
mately  1  in  1500.  The  surface  of  the  ice  shelf  be¬ 
yond  station  11  sloped  approximately  1  in  606, 
while  the  bottom  of  the  ice  shelf  had  a  slope  of 
about  1  in  137. 

If  the  brine  layer  along  the  main  profile  were  to 
maintain  a  slope  of  1  in  830  beyond  station  9  and 
could  move  instantly  inland  to  intersect  the  “im¬ 
permeable”  0.82-Mg/mJ-density  firn/ice  transi¬ 
tion  horizon,  it  would  do  so  at  a  maximum  dis¬ 
tance  of  about  14  km  inland  from  the  1977  shelf 
edge.  However,  this  is  not  likely  to  occur  under  ex¬ 
isting  infiltration  conditions  because  as  the  brine 
moves  inland  it  will  encounter  denser  and 
therefore  less  permeable  firn  as  well  as  lower  tem¬ 
peratures.  These  effects  will  slow  the  brine  infil¬ 
tration  rate  and  probably  limit  its  inland  penetra¬ 
tion  in  the  area  of  our  profile  to  less  than  1 1  km. 

A  rate  of  bottom  melting  of  the  order  of  0.9  to 
1 .2  m  of  ice  per  year  has  been  estimated  for  sta¬ 
tions  202  and  207  (see  Fig.  6)  by  Risk  and  Hoch- 
stein  (1967).  By  assuming  an  ice  bottom  ablation 
rate  of  0.9  m/yr  it  would  take  some  69  years  be¬ 
fore  the  brine-infiltrated  firn  at  station  24  be¬ 
comes  exposed  to  the  sea.  Similarly,  if  bottom 
melting  were  occurring  at  1.2  m/yr  it  would  take 
only  52  years  for  the  brine-infiltrated  firn  to  reach 
the  bottom  of  the  ice  shelf. 

Typical  ice  shelf  movement  rates  are  shown  in 
Figure  6.  These  rates  are  taken  from  Heine  (1967) 
except  for  that  of  station  307,  which  was  obtained 
from  the  1972-73  to  1978-79  austral  summer  sur¬ 
vey  data  on  file  at  the  New  Zealand  Scott  Base, 
Ross  Island,  Antarctica.  In  1978-79,  station  307 
was  moving  westward  at  84  m/yr.  The  survey  data 
also  indicate  that  the  ice  movement  at  this  station 
is  slowing  down.  If  it  is  assumed  that  the  ice  shelf 
movement  in  the  area  of  station  307  remains  con¬ 
stant  at  about  80  m/yr  during  the  next  50  to  70 
years,  and  that  the  ice  shelf  cross  section  depicted 


in  Figure  7  remains  unchanged,  then  it  is  possible 
to  make  an  approximation  of  the  general  area  of 
ice  shelf  composed  of  brine-soaked  firn  and  saline 
ice.  At  a  forward  movement  of  80  m/yr  and  a  bot¬ 
tom  ablation  rate  of  0.9  m/yr,  we  estimate  that  the 
saline  ice  near  station  307  would  reach  the  bottom 
of  the  shelf  some  5520  m  ahead  of  its  present  posi¬ 
tion  and  would  have  moved  downward  along  path 
C,  as  shown  in  Figure  7.  Similarly,  if  the  bottom 
ablation  rate  were  increased  to  1 .2  m/yr,  then  the 
brine-infiltrated  ice  would  move  forward  some 
4160  m  along  path  D  before  reaching  shelf  bot¬ 
tom.  Path  D  intersects  the  firn/ice  transition 
depth  at  a  distance  of  approximately  10,700  m  in¬ 
land  on  the  cross  section.  If  path  C  were  similarly 
extended,  it  would  intersect  the  firn/ice  transition 
at  a  distance  of  about  15  km  inland.  From  earlier 
considerations,  this  latter  distance  would  appear 
excessive,  suggesting  that  ice  shelf  bottom  melting 
is  occurring  at  a  rate  closer  to  1 .2  m/yr  rather  than 
0.9  m/yr  in  the  immediate  area  of  the  cross  sec¬ 
tion.  This  would  be  in  keeping  with  measurements 
by  Paige  (1969),  who  determined  from  ice  move¬ 
ment  records  and  drill  hole  measurements  a  shelf¬ 
thinning  rate  of  1 .06  m/yr,  which  when  combined 
with  the  annual  surface  accumulation  of  0.27 
m/yr  would  require  a  total  bottom  melting  of  1.33 
m/yr  for  an  area  located  about  0.5  km  south  of 
the  southern  end  of  our  profile  line.  This  ablation 
rate  was  subsequently  confirmed  by  Hoffman 
(1974). 

Brine  infiltration  mechanisms 
at  inland  boundary 

The  above  analysis  does  not  consider  vertical 
brine  migration,  either  by  diffusion  of  brine  along 
crystal  boundaries  (see  Eide  and  Martin  1975),  or 
by  intergranular  vein  flow  as  described  by  Nye  and 
Frank  (1973)  and  Raymond  and  Harrison  (1975). 
As  noted  earlier,  it  was  thought  that  brine  is  un¬ 
likely  to  infiltrate  much  beyond  the  firn/ice  transi¬ 
tion  because  of  lack  of  permeable  pore  space. 
However,  between  18  January  1977  and  24  No¬ 
vember  1978  the  brine  terminus  was  found  to  have 
migrated  27.5  m  closer  to  station  307,  at  an  aver¬ 
age  rate  of  0.041  m/day.  From  24  November  1978 
to  20  January  1981  the  terminus  had  moved  an  ad¬ 
ditional  30.5  m  to  a  point  18  m  beyond  station 
307,  at  an  average  rate  of  0.039  m/day.  If  this  rate 
was  maintained  then  the  brine  could  be  expected 
to  intersect  drill  hole  F  (located  37  m  beyond  sta¬ 
tion  307)  sometime  during  May  or  June  1982. 
Core  data  from  drill  station  E  revealed  that  the 
brine  was  migrating  within  ice  with  a  density  of 
0.85  Mg/m’  and  a  temperature  of  -  16 °C. 


II 


Figure  II.  Graphic  record  of  24  Novem¬ 
ber  1978. 


The  observed  brine  velocity  of  0.04  m/day  is  at 
least  three  orders  of  magnitude  greater  than  sim¬ 
ple  diffusion  theory  (through  the  ice  lattice)  would 
predict  (see,  for  example,  Hoekstra  et  al.  1965  and 
Seidensticker  1966).  Ice  with  a  density  of  0.85 
Mg/m’  would  be  permeable  to  brine  only  in  the 
presence  of  intergranular  veins  such  as  exist  in  ice 
in  temperate  glaciers.  However,  such  intergran¬ 
ular  connections  are  not  likely  to  exist  in  ice  at 
temperatures  of  -  16  °C.  One  possible  explanation 
of  continued  migration  of  the  brine  is  the  dissolu¬ 
tion  of  ice  by  the  concentrated  brine  as  it  is 
transported  into  deeper  and  warmer  parts  of  the 
ice  shelf  through  the  process  of  surface  snow  ac¬ 
cumulation  and  bottom  ice  ablation  deflecting  the 
ice  shelf  downward.  However,  considering  the 
magnitude  of  the  movement,  in  relation  to  the 
small  change  in  temperature  encountered  by  the 
brine,  it  is  unlikely  that  this  mechanism  alone 
could  account  for  the  observed  rate  of  lateral 
brine  migration.  Preliminary  field  observations  of 
thin  sections  of  brine-soaked  ice  showed  the  brine 
to  be  located  in  bubble-like  vesicles.  Unfortur~te- 
ly,  an  accidental  defrosting  of  the  McMurdo  Sta¬ 
tion  coldroom  destroyed  all  the  samples  before  de¬ 
tailed  investigations  of  the  degree  of  interconnec¬ 
tivity  of  vesicles  and  air  bubbles  in  the  brine- 
soaked  ice  could  be  completed. 

An  interesting  fact  that  may  be  related  to  a 
change  in  the  mechanism  of  brine  migration  was 
the  change  observed  in  the  attitude  of  the  brine 
layer  as  it  passed  through  the  region  of  the  firn/ice 
transition  (Fig.  11).  An  inspection  of  Figure  11 
shows  that  prior  to  reaching  the  firn/ice  transi¬ 
tion,  the  top  of  the  brine  layer  cuts  obliquely 
across  the  firn  layering.  However,  just  beyond  the 
transition  the  top  of  the  brine  layer  changes  atti¬ 
tude  to  conform  with  the  stratification  and  this 
situation  continues  to  the  brine  terminus. 

Confirmation  of  brine  depths  by  drilling 

The  calculated  depth  of  the  brine  layer,  as  deter¬ 
mined  from  the  1977  radio  echo  data,  was  verified 
by  direct  measurement  in  several  drill  holes  that 
penetrated  the  top  of  the  brine  layer.  For  example, 
using  the  January  1977  radio  echo  data  we  calcu¬ 
lated  a  depth  of  49  m  to  the  brine  layer  at  drilling 
station  E,  located  9560  m  in  from  the  shelf  edge. 
Coring  in  December  1978  revealed  a  depth  of  50.4 
m.  During  the  preceding  two  winters  about  0.9  m 
of  new  snow  had  accumulated.  Adding  this  accu¬ 
mulation  to  the  49-m  depth  gives  a  difference  be¬ 
tween  the  calculated  and  measured  brine  depths  of 
0.5  m,  or  less  than  1%.  Drilling  was  also  under¬ 
taken  in  1974  and  1977.  This  drilling  included 


three  holes  (B,  10  and  10b)  from  which  cores  were 
obtained  to  a  maximum  depth  of  22.5  m,  and  four 
other  holes  that  were  drilled  to  the  bottom  of  the 
ice  shelf  to  measure  ice  shelf  thickness.  The  four 
core  holes  drilled  in  1978  and  designated  stations 
C,  D,  E  and  F  in  Figure  7  were  originally  located 
3590,  7370,  9560  and  9700  m  respectively  in  from 
the  shelf  edge.  The  first  three  holes  (C-E)  encoun¬ 
tered  the  brine  layer  at  depths  of  19.4,  33.8  and 
50.4  ni,  respectively.  The  brine  layer  was  located 
entirely  in  permeable  firn  at  holes  B,  C,  and  D;  at 
drill  hole  E  the  brine  was  located  entirely  in  ice. 
Station  A  (Fig.  7),  located  about  5  m  from  the 
shelf  edge,  included  a  pit  excavated  2.5  m  to  sea 
level,  which  at  this  site  coincided  with  the  top  of 
the  saline  ice  of  the  brine  zone.  A  hole  was  then 
cored  approximately  6  m  to  the  bottom  of  the 
shelf  (Fig.  12).  The  core  consisted  of  brine-infil¬ 
trated  firn,  now  converted  to  saline  ice.  The  bot¬ 
tom  1-  to  2-cm  section  of  core  featured  a  dendritic 
crystal  plate  type  structure,  indicating  that  sea¬ 
water  was  freezing  onto  the  bottom  at  this  loca¬ 
tion.  The  deepest  hole,  located  at  station  F,  37  m 
beyond  station  307,  reached  a  depth  of  59  m.  Core 
data  obtained  from  the  1977  and  1978  holes  are 
listed  in  Tables  2-8. 


Figure  12.  A  ugering  into  saline  snow  ice  of  (he 
brine  zone  at  bo  t  tom  of  snow  pit  at  station  A. 
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Table  2.  Data  on  firn  from  hole  augered  on  high  side  of  brine  step  at  station  10, 
January  1977. 


Measured  Calculated 


Depth 

Temp. 

Density 

Melt 

salinity 

Brine 

salinity 

Brine 

vol. 

Air 

vol. 

Porosity 

Density * 

(m) 

(°C) 

(Mg/m') 

r/»> 

i'/m) 

r/„) 

("/«.) 

(Mg/m') 

17.15 

-15 

0.882 

15.2 

177 

59.8 

56.1 

115.9 

0.812 

17.32 

-15 

0.892 

13.0 

177 

51.7 

43.1 

94.8 

0.832 

17.47 

-  15 

— 

— 

— 

— 

— 

— 

— 

17.59 

-  14.75 

0.900 

12.2 

175 

50.1 

33.7 

83.8 

0.842 

17.82 

— 

0.894 

16.0 

175 

64.5 

44.1 

108.6 

0.819 

18.0 

— 

0.878 

14.0 

175 

63.4 

61.2 

124.6 

0.804 

18.10 

— 

— 

— 

— 

— 

— 

— 

— 

18.21 

— 

0.883 

15.5 

174 

62  5 

55.4 

117.9 

0.811 

18.35 

— 

0.896 

16.0 

174 

65.4 

42.0 

107.4 

0.820 

18.50 

0.883 

12.6 

174 

51.4 

52.4 

103.8 

0.822 

18.65 

— 

— 

— 

— 

— 

— 

— 

— 

18.73 

-  14.25 

0.887 

13.0 

172 

53.3 

48.5 

101.8 

0.823 

•  Brine-free  value. 


Table  3.  Data  on  firn  from  hole  augered  on  low  side  of  brine  step  at  station  10b, 
January  1977. 


Measured  Calculated 


Depth 

Temp. 

Density 

Melt 

salinity 

Brine 

salinity 

Brine 

vol. 

Air 

vol. 

Porosity 

Density * 

(m) 

(°C) 

(Mg/m') 

("/*) 

(’/«,) 

('/«,) 

r/x) 

c/„> 

(Mg/m') 

Table  3  (cont’d) 


Measured  Calculated 


Depth 

On) 

Temp. 

ro 

Density 
(Mg /in') 

Melt 

salinity 

C/m) 

Brine 

salinity 

C/m) 

Brine 

vol. 

C/m) 

Air 

vol. 

C/m) 

Porosity 

C/m) 

Density * 
(Mg  /in') 

19.38 

19.65 

— 

0.711 

— 

— 

— 

~ 

— 

— 

19.86 

-13.5 

0.696 

— 

— 

— 

— 

— 

— 

20.11 

-13.5 

0.694 

— 

— 

— 

— 

— 

— 

20.40 

— 

0.704 

— 

— 

— 

— 

— 

20.76 

-13.25 

0.703 

— 

— 

— 

— 

— 

— 

21.08 

-  13.25 

0.714 

— 

— 

— 

— 

— 

— 

21.45 

-  13 

0.715 

— 

— 

— 

~ 

— 

— 

21.71 

-  13 

0.727 

— 

— 

— 

— 

— 

— 

21.80 

— 

0.714 

— 

— 

— 

— 

— 

21.95 

— 

— 

— 

— 

— 

— 

— 

— 

22.10 

— 

0.800 

15.0 

162 

68.1 

144.0 

214.1 

0.732 

22.17 

-12.75 

0.874 

25.0 

162 

109.4 

75.6 

185.0 

0.747 

22.22 

— 

0.901 

30.0 

162 

135.3 

52.6 

187.9 

0.746 

22.35 

12.75 

0.903 

30.0 

162 

135.6 

50.5 

186.1 

0.748 

*  Brine-free  value. 

Table  4.  Core  data  from  station  B,  December  1978. 


Measured  Calculated 


Depth 

Temp. 

Density 

Melt 

salinity 

Brine 

salinity 

Brine 

vol. 

Air 

vol. 

Porosity 

Density • 

On) 

rev 

(Mg/m') 

C/m) 

<Vm) 

C/m) 

dm) 

dm) 

(Mg/m') 

0.05 


0.35 

— 

0.395 

— 

—  —  — 

— 

— 

0.45 

— 

0.418 

— 

—  —  — 

— 

— 

0.55 

— 

0.413 

— 

—  —  — 

— 

— 

0.65 

— 

0.395 

— 

—  —  — 

— 

— 

0.70 

— 

— 

— 

ice  lens  — 

— 

— 

0.75 

— 

0.403 

— 

—  —  — 

— 

— 

0.85 

— 

0.439 

— 

—  —  — 

— 

— 

0.90 

-  17.5 

— 

— 

—  —  — 

— 

— 

0.95 

— 

0.398 

— 

—  —  — 

— 

— 

1.15 

— 

0.402 

— 

—  —  — 

— 

— 

1.25 

— 

0.404 

— 

—  _  _ 

— 

— 

1.45 

— 

0.395 

— 

ice  lens  — 

— 

— 

1.65 

— 

0.377 

— 

—  —  — 

— 

— 

1.75 

-  17.0 

— 

— 

ice  lens  — 

— 

— 

1.85 

— 

0.392 

— 

—  —  — 

— 

— 

1.95 

— 

0.393 

— 

—  —  — 

— 

— 

2.05 

— 

— 

— 

l-cm-thick  ice  layer 

— 

— 

2.15 

— 

0.432 

— 

—  —  — 

— 

— 

2.25 

— 

0.411 

— 

—  —  — 

— 

— 

2.45 

— 

0.390 

— 

—  —  — 

— 

— 

2.55 

-  17.0 

— 

— 

—  —  — 

— 

— 

2.65 

— 

0.373 

— 

—  —  — 

— 

— 

2.75 

— 

0.353 

— 

—  —  — 

— 

— 

2.85 

— 

0.439 

— 

—  —  — 

— 

— 

2.95 

— 

0.436 

— 

—  —  — 

— 

— 

3.15 

— 

0.431 

— 

—  —  — 

— 

— 

3.25 

-  17.2 

— 

— 

—  —  — 

— 

— 

3.35 

— 

0.448 

— 

—  —  — 

— 

— 

3.45 

— 

0.466 

— 

—  —  — 

— 

— 
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Table  4  (cont’d). 

Measured 

Calculated 

Depth 

Jnt) 

Temp. 

(°C) 

Density 
(Mg /m') 

salinity 

('/«,) 

salinity 

<-/*.> 

vol. 

<y,.d 

vol. 

("/«/ 

Porosity 

(°/m) 

Density  * 
(Mg/m') 

9.70 

-3.5 

0.864 

_ 

61 

_ 

_ 

_ 

— 

9.80 

— 

0.882 

2.2 

61 

29.2 

42.9 

72.1 

0.851 

9.90 

_ 

0.882 

3.0 

61 

39.8 

44.4 

84.2 

0.840 

10.00 

_ 

0.870 

3.0 

61 

39.6 

48.8 

88.4 

0.836 

10.10 

_ 

0.871 

2.5 

61 

33.2 

55.4 

88.6 

0.836 

10.20 

_ 

0.884 

3.1 

60 

41.8 

42.5 

84.3 

0.840 

10.30 

— 

0.889 

3.1 

60 

42.0 

37.1 

79.1 

0.845 

10.40 

— 

0.889 

— 

— 

— 

— 

— 

— 

10.50 

-3.4 

0.888 

_ 

— 

— 

— 

— 

— 

*  Brine-free  value. 


Table  S.  Core  data  from  station  C,  December  1978. 


Measured 


Calculated 


Density 

salinity 

salinity 

(Mg /in') 

r/«,> 

(V„) 

Porosity 

CU) 


Density * 
(Mg /in') 


20.87  - 


•  Brine-free  value. 


ice  lenses 


brine  layer 


Table  6.  Core  data  from  station  D,  December  1978. 


Measured  Calculated 


Depth 
( m ) 

Temp. 

CC) 

Density 

(Mg/m') 

Melt 

salinity 

(°/«,) 

Brine 

salinity 

('/«,) 

Brine 

vol. 

r/ocj 

Air 

vol. 

(*/«,) 

Porosity 

("/«,) 

Density 
(Mg /in'. 

1 

0.435 

2 

— 

0.391 

— 

— 

— 

— 

— 

— 

3 

— 

0.393 

— 

— 

~ 

— 

— 

— 

4 

— 

0.436 

— 

— 

— 

— 

— 

— 

5 

— 

0.527 

— 

— 

— 

— 

— 

— 

6 

— 

0.487 

— 

— 

— 

— 

— 

— 

7 

— 

0.490 

— 

— 

— 

— 

— 

— 

8 

— 

0.555 

— 

— 

— 

— 

— 

— 

9 

— 

0.535 

— 

— 

— 

— 

— 

— 

10 

— 

0.573 

— 

— 

— 

— 

— 

— 

II 

— 

0.614 

— 

— 

— 

— 

— 

— 

12 

— 

0.571 

— 

— 

— 

— 

— 

— 

13 

-2! 

0.586 

— 

— 

— 

— 

— 

— 

14 

— 

0.607 

— 

— 

— 

— 

— 

— 

15 

— 

0.625 

— 

— 

— 

— 

— 

— 

16 

— 

0.648 

— 

— 

— 

— 

— 

— 

17 

— 

0.654 

— 

— 

— 

— 

— 

— 

18 

— 

0.671 

— 

— 

— 

— 

— 

— 

19 

— 

0.662 

— 

— 

— 

— 

— 

— 

20 

— 

0.668 

— 

— 

— 

— 

— 

— 

21 

— 

0.684 

— 

— 

— 

— 

— 

_ 

22 

-19 

0.704 

— 

— 

— 

— 

— 

— 

23 

— 

0.694 

— 

— 

— 

— 

— 

_ 

24 

— 

0.719 

— 

— 

— 

— 

— 

— 

25 

-18 

0.701 

— 

— 

— 

— 

— 

— 

26 

— 

0.718 

— 

— 

— 

— 

— 

— 

27 

— 

0.740 

— 

— 

— 

— 

— 

— 

28 

-16.5 

0.756 

— 

— 

— 

— 

— 

— 

29 

— 

0.747 

— 

— 

— 

— 

— 

— 

30 

— 

0.741 

— 

— 

— 

— 

— 

— 

31 

— 

0.742 

— 

— 

— 

— 

— 

32 

— 

0.764 

— 

— 

— 

— 

— 

— 

33 

— 

0.766 

— 

— 

— 

— 

— 

— 

33.8 

— 

0.766 

0.4 

brine  layer 

— 

— 

— 

33.92 

— 

0.781 

— 

— 

— 

— 

— 

— 

34.08 

— 

0.818 

1.8 

177 

6.6 

111.7 

118.3 

0.810 

34.26 

— 

0.824 

4.4 

177 

16.2 

107.7 

123.9 

0.805 

34.46 

-15 

0.848 

5.3 

177 

20.0 

82.6 

102.7 

0.825 

34.61 

— 

0.868 

13.1 

177 

50.7 

68.9 

119.6 

0.809 

34.77 

— 

0.892 

14.6 

177 

58.1 

44.8 

102.9 

0.824 

34.92 

— 

0.876 

8.1 

177 

31.6 

55.2 

86.8 

0.839 

35.10 

— 

0.856 

10.8 

177 

41.2 

79.5 

120.7 

0.808 

35.27 

— 

— 

13.4 

— 

— 

— 

— 

— 

35.99 

— 

— 

9.5 

— 

— 

— 

— 

36.6 

— 

— 

11.0 

— 

— 

— 

— 

— 

♦  Brine-free  value. 
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Table  7.  Core  data  from  station  E,  December  1978 


Measured  Calculated 


Melt 

Brine 

Brine 

Air 

Temp. 

Density 

salinity 

salinity 

vul. 

vol. 

Porosity 

Density  * 

rc> 

(Mg/m') 

C/oo) 

r/c) 

("/<.) 

("/«,) 

(•/»> 

(Mg/m') 

c 


Table  7  (cont’d).  Core  data  from  station  E,  December  1978, 


Measured 


Calculated 


Depth 

On) 

Temp. 

rc> 

Density 
(Mg/ in') 

Melt 
salmit , 

dm) 

Brine 

salinity 

dm) 

Brine 

vol. 

C/m) 

Air 

vol. 

C/m) 

Porosity 

dm) 

Density * 
(Mg /in') 

50.37 

brine  in  core 

50.5 

— 

0.855 

2.0 

186 

7.2 

71.9 

79.1 

0.847 

50.6 

— 

0.866 

1.3 

186 

4.7 

59.3 

64.0 

0.860 

50.7 

— 

0.869 

3.0 

186 

1 1 .0 

57.7 

68.7 

0.856 

50.8 

— 

0.867 

2.2 

186 

8.0 

59.1 

67.1 

0.857 

50.9 

— 

0.865 

1.0 

186 

3.6 

60.0 

63.6 

0.860 

51 

— 

0.855 

— 

186 

— 

— 

— 

— 

51.1 

— 

— 

10.1 

186 

— 

— 

— 

— 

51.2 

— 

— 

7.0 

186 

— 

— 

— 

— 

51.34 

— 

0.898 

9.2 

186 

34.7 

32.6 

67.3 

0.857 

51.40 

-  16.3 

0.865 

2.7 

186 

9.8 

61.7 

71.5 

0.854 

51.55 

— 

— 

2.2 

186 

— 

— 

— 

— 

51.65 

— 

0.893 

15.0 

186 

56.3 

44.0 

100.3 

0.827 

51.75 

— 

0.870 

4.1 

186 

15.0 

57.7 

72.7 

0.852 

51.88 

— 

0.867 

1.4 

186 

5.1 

58.3 

63.4 

0.861 

52.00 

— 

— 

9.6 

— 

— 

— 

— 

— 

52.11 

— 

— 

14.6 

— 

— 

— 

— 

— 

52.20 

— 

— 

0.2 

— 

— 

— 

— 

— 

52.30 

— 

— 

2.1 

— 

— 

— 

— 

— 

52.40 

— 

— 

— 

— 

— 

— 

— 

— 

52.48 

— 

— 

3.2 

— 

— 

— 

— 

— 

•  Brine-free  value. 


Table  8.  Core  data  from  station  F,  December  1978. 


Density 

(Mg/in') 


Stratigraphic 

observations 


3-  to  4-mm-thick  ice  layer 

ice  chunks 
l-cm-thick  ice  layer 


Depth  Density 
On)  (Me/m') 


Stratigraphic 

observations 


4-ram  ice  layer 


2-cm-thick  zone  of  scattered  ice  pellets 


10-cm-thick  diffuse  dirt  zone 


ice  lens 

S-cm  zone  of  ice  lenses 


5-mm-thick  ice  layer 

thin  ice  lens 
'/i-cm  ice  layer 


1 


Table  8  (cont  d). 


Density 
(Mg /m') 

Stratigraphic 

observations 

Depth 

(m) 

Density 
(Mg/ in') 

_ 

1-cm-thick  ice  layer 

40.05 

0.795 

0.726* 

40.55 

0.823 

0.721 

42.40 

0.816 

0.724 

42.60 

0.815 

— 

thin  ice  lenses  between  29.65  and  30.00  m 

43.40 

0.832 

0.737 

44.40 

0.829 

0.736 

45.60 

0.826 

0.757 

46.40 

0.831 

0.753 

47.60 

0.834 

— 

thin  ice  lens 

49.70 

0.842 

0.754 

50.20 

0.853 

0.774 

50.40 

0.846 

0.783 

51.20 

0.856 

0.773 

51.95 

0.850 

— 

2-  to  3-cm-thick  dirty  layer 

52.40 

0.856 

0.780 

53.20 

0.850 

— 

5-mm-thick  ice  layer 

54.00 

— 

— 

4-mm-thick  ice  layer 

54.95 

0.858 

— 

5-mm-thick  ice  layer 

55.40 

0.884 

0.794 

56.60 

0.874 

— 

thin  ice  lens 

57.40 

0.874 

— 

from  this  depth  on  much  of  the  core  showed  a 

58.40 

0.873 

good  banded  structure  consisting  of  fine-  and 
coarser-grained  material 

59.40 

0.865 

Stratigraphic 

observations 


3-mm-lhick  ice  layer 


l-mm-thick  ice  lens 


Density  and  temperature  profiles 

The  depth-density  profiles  for  the  firn  at  vari¬ 
ous  drill  sites  are  shown  in  Figures  13-18,  and  the 
temperature  versus  depth  profiles  obtained  at  sta¬ 
tions  B,  10b  and  D  are  shown  graphically  in  Figure 
19.  Extrapolating  the  temperature  profiles  in  Fig¬ 
ure  19  to  a  temperature  of  -  l.8°C,  the  freezing 
point  of  the  seawater,  would  yield  apparent  shelf 
thicknesses  of  23.5  m  at  station  B,  49.3  m  at  sta¬ 
tion  10b,  and  80  m  at  station  D.  The  thickness  val¬ 
ues  for  stations  B  and  10b  agree  very  well  with  the 
values  of  ice  shelf  thickness  that  can  be  obtained 
from  Figure  7.  The  value  for  station  D  is  about  9 
m  less  than  the  indicated  depth  in  Figure  7.  This 
suggests  that  the  temperature  gradient  becomes 
steeper  with  depth  at  this  location,  and  since  there 
is  no  bottom  shear  the  change  must  be  due  to  a 
change  in  density. 

Curves  obtained  by  least  squares  analysis  of  the 
depth-density  data  obtained  at  drilling  stations  B, 
C,  D,  E  and  F  are  shown  collectively  in  Figure  20. 
Each  curve  represents  the  firn  depth-density  trend 
down  to  the  top  of  the  brine  layer  except  at  station 
F,  where  no  brine  was  encountered.  These  curves, 
together  with  the  temperature  profiles  in  Figure 
19,  reveal  that  as  one  moves  further  inland  on  the 
ice  shelf  the  thermal  and  density  regimes  of  the 
firn  change.  In  short,  at  a  depth  of  10  m  at  station 


minimum  density  value — thin  slab 
missing  from  side  of  core 


B  the  snow  is  much  warmer  and  denser  than  at 
succeeding  stations  further  inland.  This  is  due  to 
more  rapid  heat  transfer  through  the  thinner  part 
of  a  wedge-shaped  ice  shelf.  Where  the  ice  shelf  is 
thin  more  heat  escapes  to  the  atmosphere  from  the 
underlying  seawater  and  larger  seasonal  tempera¬ 
ture  variations  occur  within  the  firn.  These  ef¬ 
fects,  in  conjunction  with  lower  rates  of  snow  ac¬ 
cumulation,  accelerate  the  densification  process  in 
the  firn  as  one  moves  toward  the  ice  front. 

Ice  shelf  freeboard 

From  the  data  listed  in  Table  1  we  can  calculate 
the  effective  bulk  density  of  the  ice  shelf  (<?a)  at 
station  24  from 


Ca  ('  Dt)Cw 


where  F  =  freeboard 

ew  =  density  of  seawater  ( ~  1 .0275 
Mg/m’) 

Dt  =  total  ice  shelf  thickness. 

Given  that  the  surface  elevation  is  25  m  and  the 
keel  depth  K  is  88  m,  then  DT  =  113  m,  and  from 
eq  5,  eg  is  0.800  Mg/m’. 

From  the  values  listed  in  Table  1,  the  shelf- 
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Figure  19.  Depth-temperature  pro¬ 
files  for  stations  B,  10b  and  D. 


thickness-to-freeboard  ratio  is  found  to  be 
-  4.5:1.  This  is  somewhat  higher  than  the 
—4.15:1  ratio  that  would  be  derived  from  Gow’s 
(1963)  relationship  based  on  density  data  from  the 
258-m-thick  Ross  Ice  Shelf  at  Little  America  V. 
The  higher  value  for  station  24  is  believed  due  to 
brine  loading  and  higher  firn  densities  deflecting 
the  ice  shelf  downward.  A  deflection  due  to  brine 
loading  was  quite  evident  in  our  1977  elevation 
survey,  which  showed  a  marked  change  in  slope 
near  station  10  (see  insert  A',  Fig.  7).  This  change 
in  slope  occurred  over  such  a  short  distance  that  it 
was  first  observed  visually  and  subsequently 
verified  by  the  elevation  survey  data. 

Briae  upwdliag 

Another  interesting  finding  from  the  January 


Figure  20.  Depth-density  profiles  in  brine- 
free  firn  at  the  five  stations  (B-F)  located 
along  the  main  ice  shelf  traverse  line. 


1981  radio  echo  resurvey  was  the  observation  that 
brine  had  risen  8.7  m  inside  the  0.15-m-diameter 
borehole  at  station  E  and  0.5  m  in  the  borehole  at 
station  D.  At  station  D,  the  hole  bottom  had  0.7  m 
of  hoarfrost  underlain  by  a  firm  surface  slightly 
moistened  with  brine.  At  station  E,  the  radio  echo 
data  (Fig.  21)  revealed  at  least  two  horizons  where 
the  brine  had  migrated  laterally  away  from  the 
hole.  At  the  8.7-m  upwell  level,  the  brine  had 
moved  2  to  3  m  laterally,  while  at  the  mid-upwell 
level  it  had  moved  outward  approximately  5  m.  It 
is  not  known  whether  the  8.7-m  upwelling  repre¬ 
sents  the  maximum  hydraulic  head  in  the  brine 
layer  at  station  E,  as  upwelling  may  have  been  ar¬ 
rested  by  liquid  freezing  under  the  lower  tempera¬ 
tures  at  the  higher  elevation.  Lead  line  sounding 
indicated  a  “hard  bottom”  existed  at  station  E. 
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Figure  21.  Radio  echo  profile  in  vicinity  of  drilling  station  E  show¬ 
ing  signature  of  brine  upwelling  and  lateral  migration  measured  26 
months  after  hole  was  drilled. 


Brine  chemistry 

Chemical  analysis  of  brine-infiltrated  cores 
from  several  drill  holes  along  the  main  profile 
shows  that  freeze-fractionation  of  the  seawater  as 
it  penetrates  the  ice  shelf  preferentially  precipi¬ 
tates  virtually  all  sodium  sulfate  (mirabilite)  be¬ 
fore  the  inland  boundary  of  brine  penetration  is 
reached.  For  example,  in  samples  from  near  the 
nose  of  the  4.4-m-high  brine  wave  (on  the  high 
side  of  the  step)  the  sulfate-to-sodium  ratios  are 
essentially  the  same  as  for  seawater  (0.2S),  or 
higher  (0.S9)  where  sodium  sulfate  has  precipi¬ 
tated  locally.  However,  this  ratio  decreases  4-fold 
to  0.08  immediately  beyond  the  brine  wave  (on  the 
low  side  of  the  step)  in  samples  from  an  earlier  in¬ 
filtrated  zone  across  which  the  new  brine  wave  is 
known  to  be  migrating.  Near  the  terminal  boun¬ 
dary  of  brine  infiltration,  the  sulfate-to-sodium 
ratio  is  reduced  to  0.028,  which  is  an  order  of 
magnitude  less  than  the  seawater  value.  These 
data  confirm  the  results  obtained  earlier  by  Wil¬ 
son  and  Heine  (1964).  Concomitant  removal  of 
water  by  freezing  within  the  pore  spaces  of  the  firn 
grains  has  produced  brines  with  a  salinity  approxi¬ 
mately  six  times  more  concentrated  (S  =  186V««) 
than  the  original  seawater  (5  =  33  V„).  We  present 
a  more  detailed  discussion  of  the  brine  chemistry 
and  its  gladological  implications  in  another  report 
(Cragin  et  al.  1982). 


CONCLUSIONS 

Observations  of  a  4.4-m-high  brine  step  in  the 
McMurdo  Ice  Shelf  show  that  it  has  migrated 
about  1.2  km  since  it  was  first  discovered  in  Janu¬ 
ary  1977  and  resurveyed  in  January  1981.  The 
present  wave  is  overriding  an  older  brine-soaked 
layer.  This  migration  is  proof  of  the  dynamic 
nature  of  the  step,  which  the  evidence  indicates 
represents  the  leading  edge  of  a  brine  wave  that 
originated  at  the  ice  front  during  a  recent  major 
break-out  of  the  McMurdo  Ice  Shelf.  The  inland 
boundary  of  brine  penetration,  located  8  to  10  km 
from  the  ice  front,  is  characterized  by  a  series  of 
descending  steps  that  are  also  interpreted  as  repre¬ 
senting  the  present  terminal  positions  of  several 
separate  intrusions  of  brine  similar  in  origin  to  the 
4.4-m  brine  step.  The  location  of  the  inland  boun¬ 
dary  of  brine  percolation  is  probably  controlled  to 
a  large  extent  by  the  depth  at  which  brine  encoun¬ 
ters  the  firn/ice  transition  («43-m  depth).  How¬ 
ever,  this  boundary  is  not  entirely  fixed  by  firn 
permeability  considerations  alone,  since  a  small 
but  measurable  movement  of  brine  is  still  occur¬ 
ring  in  ice  at  the  inland  boundary.  This  continued 
migration  of  brine  is  attributed  partially  to  disso¬ 
lution  of  the  surrounding  ice  by  the  concentrated 
brine  as  it  moves  into  deeper  and  warmer  parts  of 
the  ice  shelf.  Freeze- fractionation  of  the  seawater 


as  it  percolates  through  the  ice  shelf  preferentially 
precipitates  virtually  all  sodium  sulfate  by  the  time 
the  brine  reaches  the  inland  boundary.  Conco¬ 
mitant  removal  of  water  by  freezing  in  the  pore 
spaces  of  the  infiltrated  fim  produces  residual 
brines  approximately  six  times  more  concentrated 
that  the  original  seawater.  This  brine  is  retained  in 
the  ice  in  the  form  of  liquid-filled  vesicles.  We 
conclude  that  lateral  infiltration  is  the  principal 
mechanism  of  brine  soaking  of  fim  in  the  McMur- 
do  Ice  Shelf,  and  that  the  overall  process  is  dom¬ 
inated  by  wave-like  intrusions  of  seawater  trig¬ 
gered  by  periodic  break-outs  of  the  ice  front  in 
McMurdo  Sound. 
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